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G\ . Abstract 

(N 

Measurements involving the gluon spin density, AG(x,t) = G++(x,t) — G^ (x,t), 

!>■ ■ 

(^ • can play an important role in the quantitative understanding of proton structure. 

To demonstrate this, we show that the shape of the gluon asymmetry, A(x,t) = 
AG(x,t)/G(x,t), contains significant dynamical information about non-perturbative 
spin-orbit effects. It is instructive to use a separation A(x,t) = Aq(x) + e(x,t), where 
Aq(x) is an approximately scale invariant form that can be calculated within a given 
factorization prescription from the measured distributions, Aq(x, t), q(x, t) and G(x, t). 
Applying this separation with the J z = | sum rule provides a convenient way to deter- 
mine the amount of orbital angular momentum generated by mechanisms associated 
with confinement and chiral dynamics. The results are consistent with alternate non- 
perturbative approaches to the determination of orbital angular momentum in the 
proton. Our studies help to specify the accuracy that future measurements should 
achieve to constrain theoretical models for nucleon structure. 
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1 Introduction 

Several experimental programs [1, 2, 3, 4, 5] have devised strategies aimed at providing a 
significant measure of the proton's spin weighted gluon density, 

AG(x, t) = G ++ (x, t) - G + -(x, t), (1) 

where x is the Bjorken scaling variable and t = \og[a s (Ql)}/ \og[a s (Q 2 )] is the Q 2 evolution 
variable. The interest in these measurements is often framed in terms of arguments using 
the J z — \ sum rule, [6] 

Jz = \ = \ < AE(t) > + < AG(t) > +L z (t), (2) 

where 

< AE(t) > = / dxAq(x,t) and 
Jo 

<AG(t)> = i dxAG(x,t) (3) 

Jo 

are the projections of the spin carried by all quarks and the gluons on the helicity (or z)- 
axis, respectively. Also L z (t) is the net ^-component of the orbital angular momentum of the 
constituents. In this analysis, we will not attempt to identify independent flavor components 
of L z (t) within the sum rule. Arguments concerning the content and assumptions involved 
in such separations have been presented by Jaffe and Manohar [7] and by Ji. [8] Since these 
theoretical considerations do not impact the content of the analysis in sections 1 and 2, we 
will defer to the subject of separating quark and gluon orbital angular momenta until section 
3, where we discuss specific models for proton structure. 

The t-evolution of the terms on the RHS of equation (2) adds some content to this 
discussion. It was shown in [6] that the lowest order QCD evolution equation leads to 
dH/dt = in the chiral factorization prescription, reflecting the underlying helicity conser- 
vation of QCD perturbation theory. It was also shown that for a large range of boundary 
conditions, dAG/dt > 0, so that the sum rule (2) would seem to lead to the conclusion that 

lim < AG(t) >^ +oo (4) 



t— >oo 



lim L z (t) -)■ — oo. (5) 

t—toc 



The reconciliation of this result within the framework of perturbative QCD was achieved 
by Ratcliffe [9], who showed that the structure of the DGLAP evolution equations [10] is 
based upon an iteration of perturbative processes in which a change of AG = ±1 is balanced 
by a corresponding change of L z = =1=1. The chiral structure of QCD therefore implies 
that a helicity-sensitive probe of the proton's gluon density uncovers Q 2 dependence from 
Ratcliffe resolution structures that can be considered perturbatively generated P-wave color 
correlations of the composite system. Higher order evolution [10] makes small modifications 
to Ratcliffe's result without diminishing its basic content. Recent papers by Bakker, Leader 
and Trueman [11] and by Chen and Ji, [12] have clarified the understanding of angular 
momentum sum rules by dealing with some of the complications associated with the non- 
local properties of orbital angular momentum. 

Recent DIS experiments [13, 14] have significantly lowered the measurement errors of 
the quark spin contribution (AS) to equation (2). The COMPASS collaboration analysis 
quotes a result 

< AS >= 0.30 ± O.Ol(stat) ± 0.02(evol), all data (6) 

while the HERMES collaboration analysis quotes a result 

< AS >= 0.330 ± 0.025(exp) ± O.Oll(th) ± 0.028(evol), all data. (7) 

These values can be used with the J z — \ sum rule to evaluate the impact of existing and 
potential gluon asymmetry measurements. 

Recent experimental results sensitive to AG(x,t) and the gluon asymmetry, A(x, t) = 
AG(x,t)/G(x,t) have provided important new information. Although the analysis of these 
experiments is limited in sensitivity and range of < x >, the results suggest the possibility 
that at some small x and Qq, AG(x, to) < 0. Understanding the shape of AG(x,t) for the 
whole range x G (0, 1) is important, since the asymptotic result 

lim a s {t) < AG{t) > = constant (8) 

t— >oo 

follows from the full QCD evolution equations, where the constant could be positive, negative 
or zero as t — > oo. These measurements suggest that phenomeno logical conclusions described 
in equations (5) may be in error and the asymptotic color structure of the proton is quite 
different from what has been previously supposed. If the constant in equation (8) is zero, 



then < AG(t) > and L z (t) should asymptotically approach constants, while if it is negative, 
then the signs of equations (5) should be switched. Thus, experimental evidence on AG(x, to) 
for a limited range of x and to must be combined with an extrapolation in order to specify 
the nature of < AG(t) > and L z (t) at large t. The experimental results alone are not 
conclusive. The specific approach in this paper helps to illuminate these possibilities and to 
fix onto crucial experimental results. 

From the discussion above, it is instructive to write the polarized gluon asymmetry 
using the decomposition 

A(x, t) = AG(x, t)/G(x, t) = A e (x) + e(x, t), (9) 

where Aq(x) is dependent only upon x, calculable in PQCD by the using definition 

^).[(^§^)/ ( ^M)]. no) 

The numerator and denominator on the right side of equation (10) are calculable from the 
DGLAP equations and each depends strictly upon x via the respective convolutions. [10] 
The small correction, e(x, t) describes shape-dependent differences in the evolutions of G(x, t) 
and AG(x, t) at leading order (LO) in QCD perturbation theory. Differentiating A(x, t) with 
respect to t in equation (9) and using the scale- invariance of A (x), it follows that 

-[e(x,t)G(x,t)] = 0. 

The expression (9) for A(x,t) at some initial t = to leads to an equivalent decomposition for 
AG(x,t) in the form 

AG(x, t ) = A (x) ■ G(x, to) + Ag £ (x) (11) 

where the "polarized gluon excess", Ag e (x), is given by 

Ag e (x) = e(x,t)-G(x,t) (12) 

and is t— independent. This provides a boundary condition for the partial differential equa- 
tion (10) that defines A e (x) and can be used to characterize possible different shapes for 
A(x,t) in equation (9). This boundary condition for the partial differential equation (10) 
occurs at an unphysical region in that the decomposition in (11) cannot be valid when 



G(x,to) = 0. In practice, this means that there are nontrivial constraints on the magnitude 
of Ag e (x). We will discuss these constraints when we consider explicit solutions to equation 
(10). 

In this paper, we will neglect the t-dependence of AE(i,i) at LO and concentrate on 
an alternate approach to characterizing the consequences of the t-dependence of AG(x,t). 
Using the current data for AS as input, we combine equations (2), (6) and (7) to write 

L z (t)+ < AG(t) >^ -(1- < AS >) « 0.34 ± 0.02 (13) 

in a chiral factorization prescription. The quoted error is entirely due to the data uncertain- 
ties in equations (6) and (7). The asymmetry A(x,t) is parameterized as outlined in Section 
2 and the corresponding polarized gluon obtained from equation (9) is substituted into equa- 
tion (13) to determine a value for the orbital angular momentum, L z (t). Our approach to 
the study of AG(x, t) is largely complementary to the usual global analysis determination 
discussed, for example, by Hirai and Kumano [15] and others [16, 17] where the main input 
for AG(x,t) involves measurement of the scaling violations for Aq(x,t). Since it is highly 
unlikely that future experiments sensitive to AG(x, t) will determine this distribution with 
an accuracy similar to that found in the determination of Aq or to that of G(x, t) and q(x, t), 
our method takes into account the similarities between the evolution of AG(x, t) and G(x, t) 
to provide the necessary extrapolations to all x. Using knowledge of Aq(x,t) along with 
q(x, t) and G(x, t) data sensitive to values of AG(x, t) in limited regions of x and t can then 
be used efficiently to obtain valid estimates for L z (t). In particular, as long as G(x,to) is 
large enough at small scales, to, associated with confinement and chiral dynamics, L z (t) can 
also be extrapolated to provide a measure of spin-orbit dynamics at such scales. 

The possibility that there could exist a scale-invariant form for the gluon polarization 
asymmetry, A(x, t), was first considered by Einhorn [18] on the basis of numerical studies of 
the LO DGLAP equations. [10] The work of reference [19] explores Einhorn's suggestion by 
showing that a scale- invariant form for A e (x) is a feature of the Close-Sivers Bremsstrahlung 
model. Scale-invariance occurs in that approach because the perturbative diagrams which 
determine the t-evolution of the gluon distributions are the same as those that distribute the 
spin information to the gluons at LO. The papers of references [20]- [23] have significantly 
expanded the phenomenological study of the information contained in the shape of A(x, t) 



by introducing the connection to the J z = 1/2 sum rule and by considering the possibility of 
small fluctuations around the scale-invariant form described in equations (9), (11) and (12). 

The remainder of this paper is organized as follows. In Section 2 we discuss how equa- 
tion (10) leads to a non-linear equation for A € (x) based on the DGLAP equations. We 
then generate parameterizations of A e (x) in leading-order (LO) and outline the next-to- 
leading order (NLO) calculation. These parameterizations can give a meaningful description 
of A(x, t) or AG(x,t) with suitable constraints on e(x, t) or Ag t (x) respectively. When the 
constraints are satisfied, a complete description of A € (x) is obtained from boundary con- 
ditions combined with existing measurements of the singlet quark and unpolarized gluon 
densities, q(x,t), Aq(x,t) and G(x,t). Section 3 presents numerical determination of A € (x) 
within these constraints and the observation that preliminary data on AG(x, i) from several 
measurements are within these bounds. This suggests that future measurements of AG(x, t) 
or A(x,t) with improved accuracy can be used in conjunction with our analysis to determine 
a low energy non-perturbative component of orbital angular momentum. Section IV includes 
a brief discussion of the importance of understanding the constituent orbital angular mo- 
mentum in a composite system and how the J z — \ sum rule and measurements sensitive to 
A(x, t) can provide important constraints on dynamical models for proton structure. 

2 The Scale- Invariant Gluon Asymmetry Aq(x) 

2.1 Numerical approach to the gluon asymmetry 

The solution of equations based on equations (9) and (10) was proposed by Ramsey and 
Sivers [21-23]. The calculation of the asymmetry provided a method to determine AG 
without theoretical biases on its shape. This was followed by analysis of the relation between 
the Ag e parameterizations and the corresponding range of possible L z [21]. Since new data 
for the asymmetry have been made available, this method has allowed us to refine the range 
of possible < AG > and L z consistent with this data [22, 23]. The culmination of this work 
is presented in this paper. 

In this section, we create a numerical mechanism which allows phenomenological de- 
termination of the non-perturbative quantity L z from the experimental measurements of 



AG and -£r. This is done without an a priori assumption for the polarized glue. We will 
outline a numerical approach to calculate the asymmetry-!^ connection in order to provide a 
self-consistency check to this approach and provide a range of possible values for the orbital 
angular momentum. 

The perturbative component of the polarized gluon asymmetry A e Q (x) can be calcu- 
lated from a parameterization of the correction Ag e (x) by inserting equation (11) into the 
expression (10) for A e (x). Then equation (10) can then be solved using dAG/dt and dG/dt 
given by DGLAP evolution. In kinematic regions where the DGLAP evolution equations 
are valid, equation (10) allows one to generate A^(x) by 



Al(x) 



AP gq ®Ag + AP gg g (AG) 

Pgg ® q + Pgg ® G 



(14) 



where the convolution is defined as P(x) ® Q(x) = J -rP{y)Q{x/y). Then equations (10) 

and (11) can be used to generate the corresponding asymmetry using the DGLAP equations 
via the iterative equation 



Al(x) 



AP Gq ®Ag + AP GG ® [A% ■ G + Ag e 
P Gq ®q + Pgg ® G 



(15) 



Since the A e (x) term occurs on both sides of the equation, we parametrize A e Q (x) and Ag e (x) 
subject to theoretical constraints and determine the coefficients of the parameterizations that 
satisfy equation (15). For each fixed parameterization of Ag e (x), the resulting asymmetry 
A(x,t) and corresponding AG from equation (11) can then be checked for positivity with 
the corresponding unpolarized gluon at LO. The various parametrizations of Ag f (x) were 
chosen to have integrals over all x between —0.5 and 0.5. The particular parametrization 
in the second line of Table 1 was selected to change sign, consistent with instanton models. 
This possibility was discussed by the HERMES collaboration and has not been ruled out 
by data. [1, 13, 15, 29] The process to determine the NLO corrections to the asymmetry is 
similar to that of LO with the appropriate corrections to the splitting functions and evolution 
parameters. This will be addressed shortly. 

The calculation of A e Q (x) in equation (15) has been done using the CTEQ5 [24] and 
MRST [25] unpolarized parton distributions for comparison and the polarized quark distri- 
butions from reference [26]. All are evaluated at Q\ — 1 GeV 2 for the LO distributions. The 



CTEQ5 unpolarized distributions are given by 



xu v (x) = 0.9783a; a4942 (l - x) 3 - 3705 (l + 10.0012x a85n ) 

xd v (x) = 0.5959x a4942 (l - a;) 4 - 2785 (l + 8.4187x a7867 ) 

xS total (x) = 1.2214x a0877 (l - x) 7 - 7482 (l + 3.389x) 

xG{x) = 3.3862x a261 (l-a;) 3 - 4795 (l- 0.9653a;). 



(16) 



For this set, J xG(x) dx = 0.28, representing the momentum carried by the gluons. The 
MRST LO unpolarized distributions are 

xu v (x) = 0.6051x a4089 (l - x) 3 - 395 (l + 2.078^ + U.56x) 

xd v (x) = 0.0581Lr°- 2882 (l - x) 3 - 874 (l + 34.69^ + 28.96x) 

xStotai(x) = 0.2004a;- a2712 (l-a;) 3 - 808 (l + 2.283 v /a; + 20.69x) 

xG(x) = 64.57x a9171 (l-a;) 6 - 587 (l- 3.168^ + 3.25Lr). (17) 

For the MRST set, J xG(x) dx = 0.35, for comparison to the CTEQ5 set. In each case the 
total quark contribution is then given by: 

xq to tai{x) = xu v (x) + xd v (x) + xS t otai{x). (18) 

The polarized distributions given by GGR [26] in terms of the unpolarized ones are 



xAu v (x) = [1 + 0.25(1 — x) 2 J \Jx\ 1 (xu v (x) — 2xd v (x)/3) 
xAd v (x) = [1 + 0.25(1 -x) 2 /^/x]-\-xd v (x)/3) 
AS tota i(x) = (-2.36 + 2.66y/x)xS tota i(x), 



(19) 



where 



iAE(i) = xAu v (x) + xAd v (x) + xAS tota i(x). (20) 



Although, in principle, we are sensitive to the shape of the unpolarized gluon distri- 
bution at Qq, the CTEQ and MRST only differ at small x < 0.1 and the effect on our 
results is minimal. For both unpolarized CTEQ and MRST distributions used in equation 
(15), the term in the denominator vanishes at a critical Bjorken-x value of approximately 



x c ~ 0.20 — > 0.30 at LO. To avoid this numerical problem, we take the denominator and 
terms with the asymmetry to one side, giving 



Ao 



yLO 



PL ® 9 + p 99 ® G - Ap 99 ® K • G) = AP g L q u ®Aq + AP g L g ° ® (A&) . (21) 



)iO o, A „ , A DLO 



This mitigates the problem of the point discontinuity in equation (15) and makes the deter- 
mination of A e (x) straightforward. Since the asymmetry appears in both terms on the left 
hand side of the equation, a numerical trial-and-error technique must be used. 

To establish an initial reference point for the solution of equation (21), practical con- 
straints for the asymmetry at LO should include: 

• strong positivity: \A(x, to)\ < 1 for all < x < 1, and 

• endpoint values: A(0,t ) = and A(l,t ) = 1. 

Simple positivity for the gluon asymmetry at LO requires \A(x, t)\ < 1. Since the individual 
terms on the right side of equation (11) do not represent separate observables, the calculated 
Aq(x) need not satisfy these criteria, but the full calculated asymmetry, A(x,t) must. As a 
practical matter, the strong positivity constraint enforces the restriction that Ag e (x) is small 
and provides a tight restraint on the range of Ag e (x) for which stable solutions can be found. 
The constraint on endpoint values enforces the Brodsky-Farrar constituent counting rules as 
x — ¥ 1 and the requirement that AG(x = 0,t) vanishes. In all solutions, the convolutions in 
(21) are dominated in the proton by the valence up quark as x — ¥ 1. 

To investigate possible asymmetry solutions, we start with an Aq(x) parameterization 
in the form 

A e (x) = Ax a -(B- l)x fi + {B- A)x fj+1 , (22) 

which automatically satisfies the endpoint constraints provided the exponents a and /3 are 
positive. This parameterization includes substantial flexibility in adjusting shapes while 
keeping the number of free parameters to a minimum. Using the CTEQ and MRST unpo- 
larized distributions and the GGR polarized distributions given in equations (16) through 
(20) for the calculation, we find that only parameterizations of Ag e that conform to the 
condition 

| < Ag e > | < 0.25. (23) 



lead to solutions of equations (15) and (21) that satisfy the positivity constraint. A repre- 
sentative sample of the parameterizations whose asymmetries and corresponding polarized 
gluon distribution satisfy all the constraints are summarized in Table 1. These solutions give 
integrals < AG > ranging from about -0.09 to 0.59 at these small Q\ — 1 GeV 2 values. Line 
two of the table contains the only non-monotonic parameterizations of Ag e that has been 
included in this sample. It should be noted that although some of the magnitudes of the 
integrated gluon distributions, < AG >, are less than those of < Ag e >, they occur for the 
parameterizations of Ag e whose integrals are negative or zero. Since this term is added to 
the overall polarized gluon distribution, the results are consistent with the model. 

Table 1: Gluon Asymmetry Parameters at Qq — 1 GeV 2 



A<? e 


<Ag e > 


Al 


< AG> 








3a; 1 - 5 - 3x 2 - 2 + a: 3 - 2 


0.28 


a;(a;-0.25)(l-a;) 5 





1.75a; 1 - 9 - 5a; 2 - 4 + 4.25a; 3 - 4 


-0.01 


2(1 -x) 7 


0.25 


4a; 1 - 6 - 4a; 2 - 1 + x 3A 


0.51 


-2(1 -x) 7 


-0.25 


1.75a; 11 - 1.5a; 21 + 0.75a; 3 - 1 


0.18 


18x(l -x) 7 


0.25 


2a; 1 - 8 - x 2 - 5 


0.39 


-18x(l -x) 7 


-0.25 


2.25a; 17 - x 2 - 0.25a; 3 


-0.09 


-90a; 2 (l -x) 7 


-0.25 


3.5a; 1 - 3 - 4.5a; 2 - 2 + 2a; 3 - 2 


0.25 


§x(l-x) 7 


0.125 


3.75a; 1 - 4 - 3a; 1 - 6 + 0.25a; 2 - 6 


0.40 


-9x(l -x) 7 


-0.125 


3.25a; 1 - 4 - 3.75a; 2 - 2 + 1.5a; 3 - 2 


0.25 


4.5x(l - x) 7 


0.0625 


2a; - 9 - 1.5a; 12 + 0.5a; 2 - 2 


0.59 


-4.5x(l - x) 7 


-0.0625 


2.25a; L1 - 2.25a; 1 - 9 + a; 2 - 9 


0.43 



The solutions for Aq listed in Table 1 represent the best fit coefficients of equation 
(22) subject to roundoff error in solving equation (21). These parametrizations of A € are 
used to determine forms for AG(x,t) in equation (11) that can be tested experimentally in 
processes at large t for various ranges of x. The scale invariant nature of Ag e (x) implies 
that dAg e /dt = 0, so that the accuracy by which we can extract information on L z can be 



estimated. Expanding Ag e = e(x, t) ■ G(x, t) in t: 

e(x,t)-G(x,t)=e(x,t ).G(x,t ) + d{eM ;_ G{x ' t " il-f,)- "V,,;" ' (/-/,,)"+■ 



d(e(x, t)-G(x,t)) u d 2 (e(x,t)-G (x,t)) u +N2 

~dt {t ~ to)+ W 



(24) 
The second term on the right side vanishes identically and at LO and the second order term 
in equation (24) can be neglected, since second order terms in t may arise only at NLO and 
higher. Thus, at LO we can write 

AG(x, t) - AG(x, t ) « A £ (x) ■ [G(x, t) - G(x, t )\ + C(U - U,) e . (25) 

Thus, the polarized gluon distribution can be parametrized accurately at LO. All of the 
perturbatively generated orbital angular momentum can be extracted from the J z = 1/2 sum 
rule within parton distribution uncertainties to this order. At LO, AH(x) is independent of 
t, so the amount of orbital angular momentum generated by Ratcliffe resolution structures 
in the range of scales (t, to) can be given by: 

L z (t) - L z (t ) = - f'dxAl- [G(x, t) - G(x, t )] + 0(U - U,) e . (26) 

Jo 

Solving the nonlinear equation for A € (x) then leads to what we call the Ratcliffe — subtracted 
version of the J z — | sum rule, where the perturbatively generated orbital angular momen- 
tum and the perturbative evolution of AG(x,t) are subtracted. The separation allows for 
the measurement of AG(t) to be associated with a non-perturbative component of orbital 
angular momentum generated by the effects of chiral symmetry and confinement on proton 
structure. Using equations (13) and (26) we can write 

L z (t ) = 0.34 ± 0.02- < AG(t) > + < A e ■ [G(x, t) - G(x, t )] > +0(U - U,) e . (27) 

so that measurements sensitive to AG(x, t) at some large value of t can be used to estimate 
the amount of constituent orbital angular momentum present at smaller values of t in the 
perturbative region. For this study, we will choose to = (t = 0) to correspond to Q 2 = 
1 GeV 2 . This is a region where the DGLAP evolution is valid and processes evolved by 
BFKL evolution are not crucial to this study. It is also the kinematic region where chiral 
symmetry is broken. All parton distributions are evolved to this value. 

Our preliminary studies into the numerical solutions to equation (21), as indicated in 
Table 1, demonstrate that such solutions are sensitive to both the shape of Ag e (x) and its 



normalization, |< Ag € >|. This dual sensitivity arises because the region in x where the 
denominator of equation (15) is small, x G (0.24,0.30), contributes significantly to Aq(x) 
for all x < 0.30. The nonlinear nature of this equation makes the solution sensitive to the 
convolution term, APj"° ® Ag e on the right side of (21). As indicated in the discussion 
above, the requirement in (23) limits the impact of the normalization. To further study the 
numerical shape of Ag e (x), we have separately examined solutions to (21) where the support 
of Ag e (x) was restricted to the large— a; region, x > 0.30, or to the small— a; region, x < 0.20. 

These studies indicated that Ag t (x) with support in the large— a; region, the condition 
Ag e {x) <C Aq(x,t) is required uniformly in x and t for solutions to be consistent with the 
positivity condition and with the expectation that Aq(x, t)/q(x, t) — ¥ 1 as x — ¥ 1. When this 
condition is valid, the impact of Ag e (x) on Aq(x) can be ignored compared with the errors 
in Aq(x,t). Restrictions on parameterizations of Ag e (x) in the small— a; region were found 
to be much less stringent and and as long as the normalization requirement (23) was valid, 
variations in the shape in the region did not impede the ability to find stable solutions. These 
studies demonstrate that our approach of parametrizing the shape of A(x, t) can provide an 
independent, complementary approach to the global fits [15]- [17] directed at the problem 
of examining the impact of data on AG(x,t) for the determination of < AG(t) >. In 
particular, we can formulate consistent parameterizations in which the asymmetry, A(x,t), 
changes sign. 

The phenomenological application of (27) to the determination of L z (to = 0) implies 
an understanding of the errors involved in determining the components on the right hand 
side of (27). In addition to the uncertainty associated with measurements of < AS >, which 
is known explicitly, there are errors associated with 

1. the determination of < AG(t) > from experiments with limited acceptance, 

2. the numerical determination of A e (x) from solutions to (21) and 

3. the continuation to t — using < Aq[G(£) — G(0)] > within this framework. 

We will address the issues in the first two contributions in more detail in section 3. 
The errors implied by the third point can, in principle, be minimized by solving (14) at NLO 

in QCD perturbation theory instead of the LO solutions described by (21). This discussion 

follows here. 



2.2 Next-to-leading Order Asymmetry 

While we have formulated the arguments above in terms of the solutions to equation (14) 
using the LO DGLAP equations, it is appropriate at this point to insert a few comments on 
the application of solutions to the next-to-leading order set of equations 

AG NLO (x,t) = A e ( ; NLO (x,t 1 )-G NLO (x,t) + Ag e {x) 

A e,NLo ( ^ . x _ ly dAG (x,t) dG (x,t) i 

/i {x,ii) = ^ — ) \ t=tl /{ -g- ) \t=ti\ y^o) 

to the analysis of data sensitive to the polarized gluon asymmetry. Here t\ is a reference 
t at NLO. Such phenomenological analyses have traditionally been done in the framework 
of LO perturbation theory. The non-linear equation analogous to (15) that must be solved 
numerically then can be written in the form 

[A?™ ® Aq NL °(x, t) + AP G r ® [A^ NLO (x, tl ) ■ G NLO (x, t) + Ag e {x)\] 
A d i?, ti) = ± roJVLO „ ZWTTuZ In I t,nlo „ 7TWT7T7Z ,m i ( 29 ) 



\pNLO q NLO( X) t ) + pNLO GfJVLO^ f )] 

where it is specified that the convolutions on the right side are performed with the appropri- 
ate distribution functions: q NLO (x,t), G NLO (x,t), Aq NLO (x,t), AG NLO (x,t) and the NLO 
splitting functions all evaluated at t — t\. This specification acknowledges that the splitting 
functions evaluated at NLO, namely 

AP» LO (z,t) = AP$\z) + a.(t)P$\z) 
P% LO (z,t) ee P§\z)+ a ,(t)P<p(z). (30) 

involve a t dependence associated with both the MS renormalization prescription for a s (t) 

a s (t) = e- {t - h) [l + b 1 (t-t 1 ) + -- ■} (31) 

and with the factorization prescription used in the calculations. It is important to observe 
that 

AG NLO (x,t) = AG NLO (x,t 1 )+A e ' NLO (x,t 1 )-[G NLO (x,t)-G NLO (x,t 1 )]+O(U-U oo ) e , (32) 

takes a form analogous to the LO result, except that the corrections 0(t — ti) 2 to (32) are 
necessarily smaller than those found in (25). In (32), G N (x,t) evolves in t with the full 



NLO expression while evolution of t in (31) is only approximately correct at NLO. This is, 
however, an improvement on the LO result. Preliminary numerical calculations indicate that 
the same methods discussed in section 2.1 above can be used to determine fits of the form 

A^ NLO (x, h) = A(h)x a ^ - (B(h) - l)^ (tl) + (B(h) - A(h))x^ +1 , (33) 

to parametrize the shape of Aq (x,ti). 

We have not fully explored the range of boundary conditions defined by different pa- 
rameterizations of Ag € (x) for the NLO equations. For parameterizations consistent with 
(23), we have found that the differences between extrapolations for AG(x,0) generated at 
LO and NLO occur primarily for x < 0.10. It is expected, based on the above discussion, 
that errors associated with the extrapolation in t of (25) can be improved in NLO. This will 
be tested in the future. At present, however, we will address issues involving experimental 
data only at LO. 

3 Experimental Extraction of the Gluon Asymmetry 
and the Orbital Angular Momentum 

Section II describes an approach to the determination of AG(x,t) using the decomposition 
of Eq. (11) combined with experimental information on fits to q(x,t), G(x,t) and Aq(x,t) 
to solve a nonlinear partial differential equation for Aq(x). Since it uses the measured parton 
distributions in a unique way, this approach is complementary to determination of AG(x, t) 
based upon a "global analysis" as exemplified in the fits of Hirai and Kumano [15] and as 
described in other analyses. [16, 17] 

In our approach, the function Ag e (x) describes a boundary condition on the partial 
differential equation for A e (x) that represents nonperturbative dynamical constraints on the 
DGLAP evolution of gluon helicity in a proton. The limits on Ag e (x) for which stable 
solutions of Eq. (15) can be found depend sensitively upon the parameterization of the 
input parton distributions at the specific value of t at which the separation AG(x,t) = 
Aq ■ G(x,t) + Ag t (x) is defined. We have explored those limits using the specific leading 
order (LO) fits to the input distributions, Eqs. (16)-(20) at t = (Q 2 = 1 GeV 2 ). This 
tests the assumption that DGLAP evolution can be combined meaningfully with the J z — \ 



sum rule, Eq. (2), at this scale. The evolution of AG depends upon its initial shape, so its 
determination must be accompanied by assumptions of large-x behavior and interpolations 
with the small-a; shape. Our approach to the shape provides strong theoretical constraints 
and differs from global fits [15] and models such as the simpler valon model [27], which make 
different assumptions about the large- a; behavior of AG. 

For those parameterizations of Ag t {x) for which a stable solution to Eq. (15) exists at 
t — 0, we can insert the decomposition of Eq. (11) into Eq. (13) to write 

L z (0) = 0.34 ± 0.02- < A £ (x) ■ G{x, 0) > - < Ag e (x) > (34) 

and evolve the equivalent possible forms of AG(x, 0) to compare with experimental ex- 
tractions of AG(x, t) found through processes such as prompt photon and photon + jet 
measurements. The terms in Eq. (34) use the experimental averages for < AS > from 
COMPASS and HERMES, and the CTEQ or MRST parameterizations for G(x,t Q ) from 
section II. The error quoted here is due to the uncertainties in the data, as previously men- 
tioned and the small theoretical uncertainties associated with the unpolarized distributions. 
For those asymmetry parameterizations that best agree with data, we can extract the most 
likely range for the orbital angular momentum of the constituents. 

From these results, a strictly theoretical approach to finding feasible values to L z yields 
a range allowed by the "practical" constraints of 

-0.12 << AG >< 0.62 and 

-0.28 < L z (0) < 0.46, (35) 

where the theoretical and experimental uncertainties have been included in the quoted in- 
equality values. The range is not sensitive to the unpolarized input distributions and the 
corresponding theoretical uncertainties have been included in this range. The ranges of 
< AG > and < L z > appear to be numerically limited by the constraints on < Ag e > in 
equation (23). From a theoretical point of view, this results in relatively small values of AG 
and L z , but the range is still considerable. We must use the data on the asymmetry to fur- 
ther narrow the possible values of the gluon polarization and the orbital angular momentum. 
Table 3 includes a summary of present data on the asymmetry [1-5]. 

It is virtually impossible to find parameterizations that completely agree with all data. 
Two separate analyses were performed: a x 2 calculation of all parameterizations with all 



Table 2: Measurements of the gluon asymmetry 



Experiment 


Process 


Value 


stat 


sys 


•Eavg 


Comment 


COMPASS 


open charm 


-0.49 


±0.27 


±0.11 


0.11 


Large errors 


COMPASS 


High pt hadrons 


0.016 


±0.058 


±0.055 


0.09 


Q 2 > 1 GeV 2 


COMPASS 


High pt hadrons 


-0.06 


±0.31 


±0.06 


0.13 


Q 2 < 1 GeV 2 


SMC 


High pt hadrons 


-0.20 


±0.28 


±0.10 


0.07 


Q 2 > 1 GeV 2 


HERMES 


Factorization method 


0.078 


±0.034 


±0.011 


0.204 


Consistent with RHIC 


HERMES 


Approximation method 


0.071 


±0.034 


±0.010 


0.222 


Consistent with RHIC 


HERMES 


High pt hadrons 


0.41 


±0.18 


±0.03 


0.17 


Q 2 < 1 GeV 2 



existing data and a determination of which of the asymmetry parameterizations best agreed 
to within one to two a of the data set in Table 2. Six of the parameterizations in Table one, 
corresponding to lines 1, 2, 4, 6, 7 and 9 all minimized \ 2 an d were within la for at least 
five of the seven data points. All of the parameterizations that best agree with data have 
the property that Ag e < 0. A plots showing four of the parameterizations of A(x,t ) at LO 
are shown in figure 1. 

Using only these parameterizations, the range of AG and L z that result are 

-0.12 < < AG > < 0.31 and (36) 

0.03 < L 2 (0) < 0.46, (37) 

where, again, the uncertainties are included in the extreme values. These results are con- 
sistent with data and the MIT Bag model, discussed by Chen and Ji. [12] Using present 
data, the range of AG and L z have been narrowed, but clearly, accurate data over a wider 
kinematic range can more significantly constrain both the gluon polarization and the orbital 
angular momentum of the constituents. The results for the scenarios we have shown are not 
significantly different from each other. However, we do not rule out the necessity of altering 
this approach to consider more flexible parameterizations if data change significantly. We 
consider it likely that additional accurate data for the gluon asymmetry and AG over a wider 
kinematic range will allow us to further restrict our range of consistency. Two immediate 
experimental possibilities are 
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Figure 1: The existing data for the polarized gluon assymmetry with the parameterizations 
that best agree with data. 

• (1) the measurement of AG/G in a wider kinematic range of xbj for a fixed Q^ value 
and 

• (2) extraction of AG by prompt photon production or jet production in polarized pp 
collisions (RHIC). 



Our results indicate that the integrated < AG > is likely positive and small at Qq ~ 1 
GeV 2 . This is consistent with data [1, 5], chiral quark models [28] and the MIT bag model 
[12]. Although most of our parameterizations of the asymmetry are positive definite, the 
one in the second line of Table 2 changes sign and is consistent with data. This possibility 
has been discussed by others and is not ruled out by present data [12, 15, 29]. Many of our 
parameterizations give a gluon polarization consistent with zero, in agreement with much of 
the data from RHIC. [3, 4] Clearly there is still work to be done. However, we have provided 
a mechanism for calculating the gluon asymmetry that allows extraction of information on 
both AG and L z . 

In Table 3, we outline some of the (mostly quark) models that make predictions about 
the orbital angular momentum of the constituents (nearly all L q z contributions). This in- 
cludes various quark models combined with corresponding gluon scenarios. Most fall within 
the same range as our calculations, except those predicting a large orbital angular momen- 
tum, either positive or negative. However, ours is a purely phenomenological approach that 



Table 3: Comparison of Orbital Angular Momentum Models 



Model 


L z at Ql = 1 GeV 2 


Reference 


Notes 


CQSM 


0.25-0.32 


[30] 


/ig = 600 MeV 2 , no gluons 


CQSM 


0.25 


[31] 


Jg ~ 0, L q < I 


RxQM model 


0.25 


[31] 


AG < 0.10 


xQM 


0.26 


[32] 


no gluons, all L q z 


Skyrme 


0.50 


[32, 33] 


no gluons 


Casu-Seghal 


0.42 


[34] 




IK Quark model 


-0.25 


[35] 


AG = 0.59, L z w L q 


MIT Bag 


0.18 


[12, 32] 


AG: 0.2 -0.3 


RGM+OGE+vr cloud 


0.35-0.61 


[36] 


L q contribution only 


Quark model 


-0.43 


[37] 


Maximal glue - GRSV 


Quark model 


0.025 


[37] 


Normal GRSV glue 


Lattice calculation 


0.38 


[38] 


Unquenched: no glue 




Figure 2: Orbital angular momenta as a function of x for parametrizations of the polarized 
gluon asymmetry that best agree with data. 

includes gluon contributions extracted from the asymmetry and does not rely entirely on 
present fits to data. We use the data only to constrain those models which are most viable. 
Clearly, accurate data in a wider kinematic range of x and Q 2 are necessary to distinguish 
those models that are most viable and determine the total angular momentum of the con- 
stituents. 

For illustrative purposes, we have shown a plot of the x-dependent function, L z (x, t — 0) 
in figure 2 for those parametrizations of Ag e included in figure 1. Note that they have different 
behavior at small-x, further motivating the need for accurate asymmetry data, especially at 
these values of x. 

To summarize the relation between allowable values of AG and the orbital angular 
momentum of constituents, figure 3 shows the range of AG values determined by our analysis 
and the corresponding range of L z values, bounded by the vertical lines. 



4 Conclusion 



The property of color confinement implies that, at distance scales on the order of hadron sizes, 
the effective dynamic degrees of freedom for hadronic interactions are not given directly by 
perturbative interactions of quarks and gluons, but by some variety of collective excitations 
of the objects. Computer simulations of lattice regularized QCD have proven to be an 
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Figure 3: Possible range of AG values as a function of the orbital angular momentum L z 



effective tool for understanding many aspects of the crucial distance scale [41]. In addition, 
the implications of the broken chiral symmetry of QCD associated with the light quark 
masses has been incorporated into a class of effective field theories for the pseudo-Goldstone 
bosons [42] of the broken symmetry. This process has led to descriptions of nuclear structure 
such as the Georgi-Manohar chiral quark model [40] or versions of the chiral quark soliton 
model, originated by Wess and Zumino [43] in the treatment of the Skyrme solution. [33] 

Our approach to the study of the gluon asymmetry presented in sections 2 and 3 
has applied the J z — \ rule and focussed on the ability to use information on AG(x,t) 
to estimate the orbital angular momentum L z (t) at low energy scales. This is due to the 
features of proton structure that can be described in terms of the constituent quark model 
involving non-relativistic "constituent" quarks in an s— wave L = state. In this context, 
nonzero orbital angular momentum of colored constituents as inferred from equation (13) 
is associated with the internal structure of the constituent quarks. The importance of such 
virtual corrections can also be studied in the flavor asymmetry of the qq sea, [44] but the 
knowledge of L z (t) at values of Q 2 = 1 GeV 2 provides more complete information. For 
example, the orbital angular momentum of constituents in the Georgi-Manohar model has 
been studied by Song [32] and Sivers [45]. Song's results of L z = 0.30 are based upon a 
true effective field theoretical approach and does not consider the impact of gluons. The 
approach of reference [45] estimates the impact of gluonic orbital angular momentum in 
Song's approach by requiring consistency between transverse and longitudinal spin. The 
results there are L z = 0.39 ± 0.02. 

We have described a decomposition of the scale-dependence of the gluon spin asymme- 
try, AG(x,t). This approach allows a measurement sensitive to AG(x,t) to be extrapolated 
in a manner consistent with DGLAP evolution, counting rules and the measured distribu- 
tions Aq(x,t), q(x,t) and G(x,t) to help complete a more comprehensive dynamical picture 
of nucleon spin structure. The numerical work with the leading-order approximation to this 
equation indicates that the range of stability is acceptable. Furthermore, existing experimen- 
tal measurements, within reported errors suggest that AG(x, t) lies within this range. To 
illustrate the value of this decomposition, we have considered the orbital angular momentum 
inferred from the J z = | sum rule. 
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